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1.  INTRODUCTION 


If  a  gravimetric  quantity  X  is  a  function  of  £(ip.  o)  of  the  form 

X(P)=fl  ff  AW  Uv.  a)  t/o  <1.11 

o 

where  £(i^.  o)  is  some  measured  quantity  related  to  the  gravitational  potential,  then  failure  to  integrate  equa¬ 
tion  (1.1)  over  the  entire  sphere  introduces  an  error  6X  into  the  computation.  If  the  region  of  integration  is 
confined  to  a  spherical  cap  centered  on  the  computation  point  P,  the  error  is  usually  denoted  as  the  cap  trun¬ 
cation  error.  Molodenskii  et  al.  (1962)  and  recently  Dickson  (1979)  have  given  a  method  for  minimizing  d.V 
for  the  Stokes  integral.  The  method  assumes  (hat  the  first  M  harmonics  of  &f>(i p.  a)  are  known  a  priori  so  that 
one  may  replace  S(tp)  with 


Su(tfj)  =  S(ip)  -  l. 


(1.2) 


where 


M 

=  X  a"  Pm  (cos  t^).  II  J) 

*=n 

The  functions  PA  cos  i^>)  are  the  Legendre  polynomials  and  the  coefficients  a ,  are  determined  hv  minimizing 
the  integral 

J  [S(yj)  -  L,,(v>)l!  sin  V dy  (14) 

In  this  report  this  method  will  be  used  to  reduce  the  cap  truncation  error  for  the  Stokes  equation,  then 
modified  and  applied  to  the  Veiling  Meinesz  integral. 


2.  APPLICATION  OF  THE  MOLOOENSKII  PROCEDURE  TO  THE  STOKES  INTEGRAL 


Consider  first  the  Stokes  equation 

N{P)  =  ^ _ j  \  S(v)  Ag(ig,  a)  sin  d\p  do 

4nG  "  " 

having  the  kernel  function 

•S(ig)  =  X  +  *  P„(cosi*»). 

” * 1  n  -  1 


(2.1) 


(2.2) 


If  the  integration  of  equation  (2. 1 )  is  extended  only  up  to  a  spherical  distance  t p„  from  the  point  P,  then  equa¬ 
tion  (2.1)  may  be  written  as 


N(P)  =  ^  f  S(if>)  Ag(ig,  a)  do  +  AN, 


4rrG 

where  do  is  the  incremental  surface  element  and  AN,  is  the  cap  truncation  error  defined  as 

AN,  =  _ f  f  S(ip)  Ag(ig,  o)  dc. 

4nG  ° 


(2.3) 


(2.4) 


The  subscript  /  follows  the  usage  in  [Fell,  1978].  Following  the  practice  [Heiskanen  &  Moritz,  1967]  of  defin¬ 
ing  the  truncation  function  as 


S(v)~ 


0  if  0  <  ip  <  ip„ 

S(ip)  if  ifj„  <  ip  <  n 
and  expanding  it  into  a  series  of  Legendre  polynomials 


S(ip)=  X  +  1  Q.  PJcostp) 

2 

allows  the  truncation  error  to  be  expressed  as 


AN,  = 


2n  +  1 


e-XX.  Ag(ig.  a)  P„( cos  r+;)  r/o. 


4rrG  ""  2 

The  double  integral  is  equal  to  4nAg„/(2 n  +  I),  so  that 

AN,  =  JL  X  £>„  Ag„ 

2G 

where  Ag„  is  the  rt’th  degree  harmonic  from  the  expansion 

Ag  =  X  Ag„(9.  A) 


(2.5) 


(2.6) 


(2.7) 


(2.7) 


(2.8) 


and  the  Q.  arc  the  Molodcnskii  coefficients  defined  as 


2 


Q „  =  Qnityo) s  J  S  (ip)  #\(cos  41)  sin  4^/4; 

=  f  -S'( )  P„(eos  ip)  sin 

The  rms  error  due  to  neglecting  distant  /ones  is  the  square  root  of 

6N‘  =  _*1  X  g!c. 

4G2  "=2 

where  the  <•„’ s  are  the  anomaly  degree  variances  [Heiskanen  and  Moritz,  1967). 
Under  the  assumption  that  the  first  M  harmonies  of  Ag  are  known,  define 

M 

Ag  =  Ag  -  X  Ag„  =  X  Ag„ 

»r=2  «=M+I 

Then  equation  (2.1)  may  be  written  as 

M 

N  =  JL_  X  JU.  S(V)  do  +  f T  £;  S(v)  rfo 

4nG  "=2  <>  4rrG  ■> 

where  the  integration  is  carried  out  over  the  sphere. 

The  first  integral  in  equation  (2.12)  reduces  to 

M 

*  Z  =  A,  f  A',  +  ...  +  ,V„ 

G  "=2  n  -  1 

where  Ag; . Ag„  are  assumed  given. 

Let 


.S\,(4>)  =  S(v)  -  L„(tp) 


where 

M 

/  „(ud  =  a1,'  AMcos  tg) 

*=IP 


and  replace  .S'( 4/ )  by  .V„(g>)  in  the  second  integral  of  equation  (2.12). 

From  the  orthogonality  properly  of  Legendre  polynomials,  the  second  integral  is  zero  for  k 
S(\p)  can  be  suitably  modified  in  this  fashion  without  affecting  the  value  of  A.  Thus 


A  -  A,  *  ...  *  A',,  t  ^  ff  Ac  .S\,(i*>)  <7o 
4nG  ' 

a,  +...♦  a,, f  ~%r 
4nO 

where  d.V„  is  the  cap  truncation  error  under  the  assumptions  made  above. 


d.A„  =  R  J  J  Ag  S„M  r/o. 
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(2.9) 

(2.10) 

(2.11) 

(2.12) 

(2.13) 

(2.14) 

(2.15) 

A/.  Thus 

(2.16) 


R 

4nG 


(2.17) 


Let 


= 


0  iT  0  <  tp  <  ip„ 

S„( ip)  it'  ip„  <  ip  <  n 


and  expand  it  into  a  series  of  Legendre  polynomials 
.S„(ip)  =  1  2^_JL!  r/I'P.fcosip) 

r-o  2 

where 

q"„  =  f/"(tp,,)  =  J  S  „(ip)  P„( cos  ip  I  sin  ipr/ip 
=  I  S\,(ip)  P„(cos  tp)  sin  tprftp. 
Using  equations  (2.14)  and  (2.15)  and  letting 
/?“*  =  I  P„(cos  ip)  P4(cos  ip)  sin  ipz/ip 
one  can  rewrite  equation  (2.20)  as 

\t 

=  Q„  -  2.  alR'l, 

k  *0 

and  the  cap  truncation  error  as 


d,V„  = 


R  ...  X  2"  "  1  </„’  /  /  Ac  P„ 
nr;  -»»  ?  "  " 


■inCl 


&f!  P„(cos  ip)  sin  ipr/ipr/o. 


(2. IX) 


(2.19) 


(2.20) 


(2.21) 


(2.22) 


(2.23) 


Since  the  integral  in  equation  (2.23)  is  zero  for  n  <  M  and  equal  to  4nA.g„/(2/t  +  I)  for  n  ^  M  t  I, 

d.V„  R  X  </.:  Ag„.  (2.24) 

20’ 

The  rms  error  due  to  neglecting  distant  /ones  is  the  square  root  of 


d.v; 


4(7 


;  1  h:fpl 


U/.' )  c. 


(2.25) 


To  evaluate  equation  (2.25)  the  coefficients  </"  must  he  computed.  This  is  done  by  applying  Schwarz's  in¬ 
equality  to  equation  (2.17), 


| JJ  Ac  S„(ip)  r/oj  <  JJ  (Ae 1 1  do  •  JJ  [.S„(ip)I /  do. 


(2.26) 


4 


and  minimizing  the  second  integral  on  the  right.  This  leads  to  the  following  set  of  M  +  I  equations  with  a 
many  unknowns 


I 


a"  P,(cos  ip) 


=  0 


for  n  =  0,  1 . M 


giving 


f”  S{\p)  P„(ip)  sin  xpdxp  = 


P„(  cos  tp)  Pi  (cos  tp)  sin  ipdip. 


Using  equation  (2.9) 

M 

Q„=  X  a'l  R'l„  for  n  =  0,  l, ....  M. 
*=0 


(2.27) 


(2.28) 


Expressing  equation  (2.28)  in  matrix  notation 

{Q)  =  [R°]  {a0} 

gives  the  solution  as 

{«"}  =  W0)"  \Q). 

Note  that  a"„  =  a"(M,  i p„),  i.e.,  a»  is  a  function  of  M,  the  number  of  lower  degree  harmonics  assumed  known, 
and  of  the  cap  size,  tp„.  Furthermore,  q"„  =  0  for  n  <  M. 


3.  APPLICATION  OF  THE  MOI.ODENSKII  PROCEDURE  TO  THE  VENING  MEINESZ  INTEGRAL 


The  two  components  of  the  deflection  of  the  vertical  are  the  north-south  component  i  and  the  east-west 
component  rg  The  Vening  Meinesz  integral 

{  U  =  —i—  X  X  S  (v)  a)i  sin  »  }  do  <-v  1  > 

4nO 

where 

S'(ig)=  —  S(tg)  =  -  X  2"  +  1  P.!(cosig)  (3.2) 

r/tg  '"•*  //  -  I 

permits  the  computation  of  the  deflection  of  the  vertical  from  gravity  anomaly.  If  the  integration  is  carried 
out  over  a  spherical  cap  of  radius  ig,„  then 

f  ".S  (tg)Ag{a"u}r/o  +  {f  }  (3.3) 

J  ^  Sin  O  1  1  (1 Y]  1  I 

I)  -'ll 


where 

<£},=  0A) 

4nG 

is  the  cap  truncation  error. 

Following  [Hagiwara,  1973]  introduce  the  function 


4M  = 


0 

\  s  (v) 


if  0  <  tg  <  Ign 
if  tg„  <  tg  <  n 


and  expand  it  into  a  series  of  associated  Legendre  functions 

4(ig)=  5-  .  r/„(tg„)  P,!(cox  tg) 

2 n(n  4  1) 

where  </„( ig„)  are  the  Cooke  coefficients  defined  by 

qA tgu)  =  !,  X  S  (>g)  f>-  (cos  tg)  sm  i gr/tg 

Furthermore,  Hagiwara  shows  that  r/„( tg„)  ;tnd  (L,(tg„)  ;,rc  related  by 

r/„(i g„)  =  -  *  **  £L(tgi>)  "  l  .V(ig„)  /’,!(cos  ig„ )  sin  tg„. 

2 

Using  equations  (3.5)  and  (3.6)  and  defining 


(3.5) 


(3.6) 


(3.7) 


(3.8) 
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(3.9) 


Q*  =  Q*(ip „)  =  - ? - q„{\p  o) 

n(w  +  1) 

=  <?.,  +  _ ! _ S(ipo)  P,',(cos  ip„ )  sin  i p„ 

n(n  +  1) 

the  cap  truncation  error,  equation  (3.4),  becomes 


|<#l  = 

l  <>i  1 1 


2n  G 


i  2”  ii-  q:Jh  fo  a«  pu cos  V)  { r;,:: }  do. 


(3.10) 


It  can  be  shown  [Hagiwara,  1973),  [dcWitte,  1967)  that  the  double  integral  in  equation  (3.10)  through  a 
change  to  latitude  and  longitude  variables  is  equal  to 


f„  f„  Ag  Pi  (cos  W)  {  ^'ns  l}  da 


4n 


2  n  +  1 


a 

Yt> 

t  a 

cos  $  ax 


Ac. 

Ac. 


(3.11) 


So  the  cap  truncation  error  becomes 


/ « i  =  _  i  y  Q* 

Ont,  w  £  «■ 


I 


a* 

a 


Ac. 


Ac. 


(3.12) 


cos  *  a\ 

and  the  rms  error  due  to  neglecting  distant  zones  is  the  square  root  of 


d0,2=£{d42  +  dr,2}  =  _J _  X  n(n  +  t)<2„*2c-„  (3.13) 

4  G‘  "=2 

where  £  is  the  expected  value  operator. 

Following  a  procedure  similiar  to  that  used  for  the  Stokes  integral,  assume  that  the  first  M  harmonics  of 
Ag  are  known. 

Let 


i„(4j)  -  .S'  (i p)  -  Ltiv) 


(3.14) 


where 

M 

Ll(ip)=  X  al  Pi  (cos  4j), 

fc=l 

and  in  equation  (3.1)  replace  S  (ip)  by  £„( tp).  This  yields 


{  4  }  =  _ 

4nO 


4n(i 


M  > 

I,  I,  A*"  s(ip)  { <s;;;; }  do 

-  f  f"a  { v.n  r }  + 1  aU  „  • 

r  "'ll  -^ti 


(3.15) 


(3.16) 
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Using  equation  (3.1 1)  the  first  integral  in  equation  (3.16)  can  be  evaluated  giving 

a 


1 


1  I 

C,  (/;  -  1) 


i 

CON  ' 


a* 

a 

a* 


A*. 

A*. 


The  last  term  in  equation  (3.16)  is  the  cap  truncation  error 

<S1„- 


Let 


Mv)  : 


4  nG 

0  if  0  <  \p  < 

'  i«(v)  if  4-„  <  ^  <  rr 


and  expand  it  in  a  series  of  associated  Legendre  functions  of  order  one, 

4  =  5!  2/1 .  +-1 —  PUcos  tp ) 

"=l  2«(«  +  1) 


where 


<7„  s  s  J  4  „(t fj)  P'„ (cos  v)  sin 

=  ^  4«M  p!,  (cos  4>)  sin  n>d\p. 

Using  equations  (3.9)  and  (3.15)  and  defining 
R  !,*  2  J  P'„  (cos  ip )  Pi  (cos  tfi )  sin  \yd\p 


ql  =  q'Awn)  =  Q*  + 


1 


I  a,’*: 


«(«  +1)  "=l 
the  expressions  for  q„(q>»)  and  £„(ip)  become 


d-  =  </., 


!,  X  a  l  R  '  „» 


_ «(«  +■  1) 


n(n  +  1) 


C??  + _ ' _  ^  a',  Rl 

n(n  +1)  **' 


(/■. 


(3.17) 

(3.18) 

(3.19) 

(3.20) 

(3.21) 

(3.22) 

(3.23) 


(3.24) 


8 


and 


£  (i p)  =  ~  </ '  p.!  (cos  ^ ). 

"='  4 

From  equation  (3.9)  we  note  that  the  </.!  are  related  to  the  Molodenskii  coefficients 


q\  ~  Q»  + _ ! _ S(ip„)  P.',  (cos  tg„)  sin  tg„  +  1  X  r/t'  /f  1» 

n(n  +  I)  n(n  +1)  *=l 


The  cap  truncation  error,  equation  (3.18)  ean  now  be  expressed  as 

IS),,--  '  t 


=  -  J_  I  r?.:  <’♦ 

2G  l  <1  At 

e»'+  ,n 

and  the  total  rms  error  due  to  neglecting  distant  zones  as  the  square  root  of 


n(n  +  l)lq.',]‘c„. 


(.3.26) 


Before  computing  the  coefficients  q'„ ,  we  need  to  determine  the  coefficients  al .  Applying  to  equation  (3.18) 
the  same  procedure  as  was  used  for  Stokes  integral  leads  to  the  condition 


— —  If  f  S  (ig)  -  5!  a  l  P'k  (cos  \p )  rfoi  =  0,  for  n  =  1, 2 . . 

)</'„  l”  **'  J  J 


which  results  in  the  linear  equations 


X"  M  „ 

-S'  (ig)  PUcos  \p )  sin  tgr/tg  =  I  [  P„'(cos  y)  PU cos  iq)  sin  i pdw,  for  n  =  1,2 . M  (3.30) 

^  U'«, 


I-  -  1,  X  a,'/i  1,.  for  n  =  1,2 . A/. 


(Expressing  equation  (3.31)  in  matrix  notation  gives 


9 


{q)  =  \  \r']  (fl'i 
giving  the  solution 

{a'}  =  2l  R'y'iq).  (3.32) 

Note  that  a'k  =  a'k(M,  v>0)  and  from  equation  (3.24) 

<7„  =  q'„  =  0  for  n  «  M. 
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4.  COMPUTATIONAL  APPROACH  FOR  KVALUATING  THF  TRUNCATION  F.RROR, 

RESULTS,  AND  CONCLUSIONS 


The  truncation  error  for  the  Stokes  and  Vening  Meinesz  integrals  due  to  neglecting  distant  /ones  was 
given  by  equations  (2.10)  and  (3.13)  respectively.  Introducing  the  modification  to  the  kernel  recommended  by 
Molodenskii  (1962)  leads  to  equation  (2.25)  for  geoid  height.  For  vertical  deflection  an  analogous  procedure 
leads  to  equation  (3.28). 

Other  methods  for  minimizing  the  truncation  error  have  been  studied  [Wong  (1969),  Meissl  (1971),  Fell 
(1978)].  In  addition  to  the  Molodenskii  procedure,  results  in  this  report  include  the  procedure  recommended 
by  Fell  (1978)  for  the  Stokes  kernel,  consisting  of  harmonic  removal  from  the  gravity  anomaly  field  alone. 
Results  from  this  procedure  are  presented  for  both  the  Stokes  and  Vening  Meinesz  integrals. 

For  this  latter  method  the  expected  truncation  error  for  geoid  height  is  given  by  the  equation 


6N‘  = 


R 2 

4 G7 


I 


0U\ 


(4.1) 


..here  the  Q„  are  the  Molodenskii  coefficients  given  by  equation  (2.9).  For  vertical  deflection  the  expected 
truncation  error  for  (his  method  is 


60',=  — ! -  X  n(n  +  1)  £>.?'<'„  (4.2) 

4G 1 

where  the  Q*  coefficients  are  given  by  equation  (3.9)  as  developed  by  Hagiwara  (1973). 

In  addition  to  equations  (4.1)  and  (4.2)  the  equations  used  to  evaluate  the  truncation  errors  are  sum¬ 
marized  as  follows: 


6N,'=  ZL_  I  tfr. 

4rrG2  »=> 

6*:,=  t  (r/lV  c„ 

4nG 1 

6d‘=  _ 1  Z  n(n  +  1)Q!'c„ 

4G  "=J 


60'  = 


1 

4G1 


51  nln 


n(n  +  l)(r/„)  c„ 


where 

Q-  =  f  ,  S(x)  P„(.\)  d\ 

M 

O  "V  II  r.o 

7"  ~  C?"  “  ft  n* 

k- 1> 

R'L  =  /\(.v)/\(v)r/.v 

Q*  +  1  .S(.v„)|/V1(.v„)  -.v„P.(.v„)l 

</l  =  Q*,  *  .....!  5  al  R I. 

n(n  *  I) 


(4.3) 

(4.4) 

(4.5) 

(4.6) 

(4.7) 

(4.8) 

(4.9) 

(4.10) 

(4.11) 


II 


PL  =  J  P'Ax)  P'Ax)  dx 

(4.12) 

S(x0)  =  1  -  6  sin  V?  +  l  -  cos  t^0  [5  -  3  ln(sin  V°  +  sin 2  V«)] 

sin  Vo  2  2  2 

(4.13) 

T 

Xo  =  COS  tp„ 

(4.14) 

All  computations  were  done  on  the  CDC-6700  computer  using  double  precision  arithmetic.  The  anomaly 
degree  variances  were  computed  using  Kaula’s  rule 

C„=  192  for  n  ^  3  (4  |5) 

n  +  1.5 

with  c2  taken  as  10  mgals2. 

The  values  of  several  constants  were  taken  as 

G  =9.798  x  io’  mgals 

R1 /Gl  =  42.3  mVmgaP.  (4.16) 

The  Molodenskii  coefficients,  £>„,  were  computed  using  the  recurrence  relations  derived  by  Paul  (1973).  To 
compute  certain  other  quantities,  the  following  identities  were  used  (See  Appendix  C  for  derivation). 

/  \P„(x)Ydx=  2n  ~  1  J  Pl,{x)dx+  U/^Cy)  +  Pl,(x)-2PAx)P..,(x)  1*  for#|>, 

2n  +  1  “  L  (2n  +  I)  Jo  (4.17) 

rh 

J  PAx)  P„(x)  dx  =  (n-k)x r  PAx)  PAx)  -  nP„-Ax)  PAx)  +  kPAx)  P,-,(x)  ]" 

(n  -  k)(n  +  k  +  1)  \u 

for  A:  and  «>1,A:>I  (4.18) 

J  \P'Ax)Y  dx  =  n(n  +  \)  J  \PAx)Y  dx+  n\PAx)  P„.Ax)  -  xP‘Ax))o  for  n  »  1  (4.19) 

J  P'Ax)  P'Ax)  dx  -  -  nk  \(n  -  k)xPAx)  PAx)  -  (n  +  I)  P.(.v)  P,.,(x)  +  {k  r  1)  P„  Ax)  PAx)  ]" 

L  (n  -  k)(n  +  k  +  I)  J 

for  n  *  k  and  n>2,  k>  2.  (4.20) 
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GEOID  HEIGHT  RESULTS 


Figure  1  presents  a  plot  of  the  classical  geoid  height  truncation  error  as  a  function  of  cap  si/e  based  on 
equation  (2.10).  This  graph  shows  that  the  expected  error  in  geoid  height  does  not  decrease  asymptotically  as 
the  area  of  integration  is  increased  as  might  be  expected  intuitively.  The  error  function  has  two  local  minima 
occurring  at  the  zeros  of  the  Stokes  kernel.  Even  if  the  cap  si/e  is  extended  to  sixty  degrees,  the  truncation  et  - 
ror  still  exceeds  10  meters. 


Figure  I.  Truncation  Error  for  Geoid  Height 


To  reduce  this  error,  two  procedures  were  introduced  with  the  resultant  truncation  errors  [Fell, 
1978]  and  <5iVw  [Molodcnskii,  1962].  These  procedures  were  evaluated  for  several  values  of  M  representing  the 
number  of  harmonics  of  the  gravity  anomaly  field  which  are  assumed  to  be  known.  Figure  2  gives  the  resul¬ 
tant  truncation  errors  for  these  procedures  for  the  case  where  M  is  6.Jn  this  case  the  truncation  error  was 
substantially  reduced  by  the  removal  of  6  harmonics  from  the  field  (dW ,„)  and  dramatically  reduced  even 
further  by  the  Molodcnskii  technique  (d/V>f).  As  an  example,  for  a  cap  si/c  of  twenty  degrees  the  classical 
truncation  error  of  14.8  meters  was  reduced  to  2.9  meters  by  the  removal  of  6  harmonics  from  the  gravity 
anomaly  field  and  reduced  to  0.5  meters  using  the  Molodcnskii  procedure  which  modifies  the  Stokes  kernel 
according  to  equation  (2.14).  Figure  3  gives  the  modified  and  original  Stokes  kernel  for  this  example,  and 
Figure  4  presents  the  truncation  errors  dA/(()  and  (W,t  as  a  function  of  M  for  a  twenty  degree  cap. 
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CAP  SIZE  it>o  (degrees) 


Figure  2.  Expected  Truncation  Error  for  Geoid  Height  Using  Method  III  and  Molodenskii’s  Procedure  {M  =  6). 


(degrees) 


Figure  3.  Stokes  Function  and  Modified  Kernel  According  to  Molodcnskii  {M  =  6,  v«  ~  20"). 


1 

1 


Appendix  A  contains  additional  information  and  results  obtained  using  these  two  procedures  for  reduc¬ 
ing  the  truncation  error.  The  polynomial  coefficients  a>(M,  )  for  the  kernel  modifications  using  Moloden- 

skii’s  procedure  arc  given  for  several  values  of  M.  Examples  of  the  modified  Stokes  kernel  .S’„(ig)  are  also 
plotted  in  this  appendix.  An  examination  of  these  figures  shows  that  for  increased  values  of  \f  or  for  larger 
cap  sizes  the  modified  kernel  according  to  Molodenskii  approaches  zero  outside  the  cap  as  expected  since  this 
was  the  aim  of  the  minimization  criterion.  Inside  the  cap  the  modified  kernel  appears  somewhat  like  the  result 
expected  of  a  cosine  taper  applied  to  S(\p).  Figures  A. 3  give  additional  examples  of  the  truncation  error  as  a 
function  of  cap  size  for  M  =  2,  4,  12  and  Figures  A. 4  give  the  truncation  errors  AN m  and  AN  as  a  function 
of  M  for  cap  sizes  of  15  and  30  degrees.  Finally  Tables  A. 2  give  the  truncation  errors  for  M  =  2.  4.  6,  8,  10 
and  12  as  a  function  of  cap  size  ranging  to  100  degrees.  These  results  demonstrate  for  most  cases  that  the 
Molodenskii  technique  is  an  extremely  effective  procedure  for  reducing  the  geoid  height  truncation  error.  The 
procedure  examined  by  Fell  (1978)  consisting  of  harmonic  removal  from  the  anomaly  data  alone  also 
significantly  reduces  the  error  of  truncation  but  not  nearly  as  well  as  the  optimized  procedure  of  modifying 
the  kernel.  The  subtraction  of  the  first  M  harmonics  of  the  Stokes  kernel  with  coefficients  (2n  +  l)/(«  -  I) 
from  .S'fgd  is  another  (not  optimized)  procedure  examined  by  Fell  (1978)  which  cannot  be  expected  to  yield  as 
good  of  results  as  with  the  Molodenskii  technique.  However  the  Molodenskii  procedure  is  weakest  when  ap¬ 
plied  with  small  M  and  small  cap  size  as  evidenced  in  Tables  A. 2.  In  these  cases  the  minimization  procedure 


i 
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VERTICAL  DEFLECTION  RESULTS 


The  truncation  error  for  vcriical  deflection  based  on  equation  (3.13)  is  given  in  f  igure  5  as  a  function  of 
cap  si/e.  F  or  instance,  using  observations  within  a  forty  degree  cap  would  result  in  an  rms  truncation  error  of 
approximately  !  second  of  arc. 


F  igure  5.  Truncation  Error  for  Vertical  Deflection 


To  reduce  this  error  two  procedures  arc  examined.  The  first  is  based  on  removal  of  harmonics  from  the 
gravity  anomaly  field  adopting  the  Vening  Mcinesz  kernel,  while  the  second  additionally  modifies  the  kernel 
according  to  equation  (3.14).  In  this  latter  approach  the  Molodcnskii  technique  for  geoid  height  is  modified 
by  changing  the  form  of  the  approximating  polynomial.  Figures  6  and  7  give  examples  of  this  modified  kernel 
for  M  equal  to  6  for  two  cap  si/cs,  20  and  40  degrees.  For  the  latter  case  the  modified  kernel  resembles  a 
cosine  taper  applied  to  the  Vetting  Meines/  kernel.  In  Figure  H  the  expected  truncation  error  for  vertical 
deflection  as  a  function  of  cap  si/e  is  illustrated  for  the  three  cases  mentioned.  Substantial  reduction  in  the 
truncation  error  is  apparent  using  Method  III  and  the  Molodcnskii  technique,  with  the  latter  producing 
better  results  for  caps  si/es  exceeding  20  degrees.  In  this  example  the  Molodcnskii  procedure  would  indicate  that  a 
cap  si/e  of  between  20  and  30  degrees  would  be  adequate  for  a  deflection  accuracy  of  0.2  seconds  of  arc.  This  in¬ 
formation  is  thus  valuable  tor  survey  design,  l  or  a  desired  accuracy  and  for  an  assumed  value  for  M  a  minimum 
cap  si/c  may  he  determined  from  graphs  analogous  to  Figure  N.  If  deflections  are  desired  in  a  limited  geographical 
region,  then  the  extent  of  the  survey  to  support  deflection  computations  in  that  region  may  be  defined. 
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Figure  6.  Vening  Meinesz  Kernel  and  Modified  Kernel  Based  on  a  Molodenskii  Approach  (M  =  6,  4/,,  = 
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4/  (degrees) 

Figure  7.  Vening  Meines/  Kernel  and  Modified  Kernel  Based  on  a  Molodenskii  Approach  ( A#  6, 


II  20  40  60  SO  100  120  140 


CAP  SIZE  <4>0  (degrees) 

Figure  8.  Expected  Truncation  Error  for  Vertical  Deflection  as  a  Function  of  Cap  Size  (M  -  6). 


Figure  9  gives  an  example  of  the  rms  deflection  error  for  a  30  degree  cap  as  a  function  of  the  number  of 
harmonics  assumed  known.  Again  for  small  caps  and  small  M  the  approximation  of  the  Vening  Meinesz 
kernel  with  the  low  degree  polynomial  will  not  be  as  effective  a  procedure,  since  the  modified  kernel  may  be 
significantly  different  from  zero  exterior  to  the  cap. 

In  Appendix  B  more  detailed  results  are  presented  for  these  two  procedures  for  minimizing  the  trunca¬ 
tion  error.  Table  B.l  gives  the  coefficients  al(M,  tp„)  for  the  modified  Vening  Meincsz  kernel  based  on  the 
Molodenskii  procedure  for  M  equal  to  2,  4,  6,  8  and  10.  Figures  B.  1  and  B.2  give  examples  ot  these  modified 
kernels  for  cap  sizes  of  20  and  40  degrees  respectively.  Figure  B.4  gives  the  expected  truncation  error  for  ver¬ 
tical  deflection  using  Method  III  and  the  Molodenskii  approach  as  a  function  of  the  number  of  harmonics 
removed  for  the  20  and  40  degree  spherical  caps.  And  finally,  in  Table  B.2  the  truncation  errors  69  r  69  and 
69  M  arc  given  for  cap  sizes  up  to  100  degrees  for  AT  values  of  2,  4,  6,  8,  10,  and  12.  Again  for  small  values  of 
M  and  these  tables  show  that  the  Molodenskii  approach  and  Method  ill  give  comparable  results.  As  the 
number  of  harmonics  removed  and  the  cap  size  increases  the  Molodenskii  approach  becomes  superior. 
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Number  of  Harmonics  kemoved  (M) 

Figure  9.  Expected  Truncation  Error  for  Vertical  Deflection  as  a  Function  of  the  Number  of  known  Harmonics  30 c) 
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This  report  has  addressed  some  relatively  simple  approaches  for  minimizing  truncation  error  for  geotd 
height  and  deflection  of  the  vertical  computations  due  to  the  application  of  the  Stokes  and  Vcning  Meinesz  in¬ 
tegrals  on  limited  geographic  areas.  The  results  clearly  demonstrate  that  major  reductions  in  these  truncation 
errors  are  possible  under  the  assumption  that  the  lower  degree  harmonics  of  the  disturbing  potential  are 
known.  The  effects  of  errors  in  these  harmonics  and  errors  in  gravity  anomaly  data  were  not  considered. 
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APPENDIX  A 

MODIFIED  STOKES  KERNELS  AND  TRUNCATION  ERRORS 


Table  A. I  Polynomial  Coefficients  a "(\f,  ip„ f  for  Stokes  Kernel  Modifications  According  to  Molodenskii 


M  =  2 


10.0 

-.18183 

-.54283 

4.10415 

15.0 

-.26676 

-.79184 

3.70801 

20.0 

-.34709 

1 .02267 

3.35591 

25.0 

-  42365 

-1.23737 

3.04532 

30.0 

-.49728 

-1.43820 

2.77224 

35.0 

-.56865 

-1.62715 

2.53235 

40.0 

-.63827 

-1.80583 

2.32154 

45.0 

-.70647 

-1.97549 

2.13608 

50.0 

-.77350 

-2.13706 

1 .97267 

55.0 

-.83947 

-2.29123 

1.82843 

60.0 

-.90444 

-2.43849 

1 .70088 

M  =  4 

if»n 

al 

</'! 

a') 

a'\ 

a" 

10.0 

-.14969 

-.44711 

4.26130 

2.47930 

1.71041 

15.0 

-.20709 

-.61574 

3.99193 

2.12596 

1.29522 

20.0 

-.25791 

-.76264 

3.76498 

1.84379 

.98873 

25.0 

-.30444 

-.89478 

3.56835 

1.61 391 

.76178 

30.0 

-.34805 

-1.01639 

3.39429 

1 .42333 

.59270 

35.0 

-.38956 

-1.13005 

3.23785 

1.26312 

.46579 

40.0 

-.42947 

-1.23736 

3.09575 

1.12696 

.36975 

45.0 

-.46807 

-1.33933 

2.96575 

1.01028 

.29645 

50.0 

-.50555 

-1.43662 

2.84625 

.90968 

.24004 

55.0 

-.54200 

-1.52965 

2.7.3605 

.82250 

.19626 

60.0 

-.57749 

-1.61872 

2.63423 

.74668 

.16200 

M  -=  6 

ah 

a", 

a" 

a" 

tf” 

a'!. 

10.0 

-.13010 

-.38876 

4.35713 

2.61052 

1.87419 

1 .4007 1 

1 .04240 

15.0 

-.17431 

-.51892 

4.14827 

2.33460 

1.54658 

1 .03965 

.66796 

20.0 

-.21263 

-.630.34 

3.97408 

2.11371 

1.29918 

.78807 

.43449 

25.0 

-.24742 

-.73017 

3.82216 

1.92909 

1.10486 

.60724 

.28720 

30.0 

-.27985 

-.82203 

3.68603 

1.77055 

.94808 

.47404 

.19298 

35.0 

-.31056 

-.90792 

3.56205 

1 .63206 

.81930 

.37409 

.13182 

40.0 

-.3.3990 

-.98899 

3.44799 

1.50977 

.71223 

.29802 

.09152 

45.0 

-.36812 

-1.06598 

3.34238 

1 .40099 

.62243 

.23950 

.06456 

50.0 

-.39533 

- 1 . 1 3936 

.3.24419 

1.30375 

.54664 

.19405 

.04627 

55.0 

-.42164 

-1.20946 

3.15264 

1.21652 

.482.36 

.  1 5850 

.03367 

60.0 

-.44709 

-1.27649 

3.06716 

1 . 1 3808 

.42761 

. 1 3049 

.02488 

Table  A.l  (Continued) 


M  =  8 


tp,i 

a‘l 

a" 

a" 

a " 

u" 

a\\ 

</*' 

i i'l 

10.0 

-.11624 

-.34748 

4.42497 

2.70348 

1.99031 

1.53757 

1.19716 

.92073 

.68826 

15.0 

-.15273 

-.45515 

4.25146 

2.47272 

1.71368 

1.22864 

.87095 

.59298 

.37456 

20.0 

-.18421 

-.54707 

4.10635 

2.28578 

1 .49933 

1 .00306 

.65064 

.39329 

.20838 

25.0 

-.21271 

-.62951 

3.97884 

2.12664 

1.32482 

.83013 

.49480 

.26672 

.11861 

30.0 

-.23922 

-.70539 

3.86380 

1.98747 

1.17874 

.69368 

.38122 

.  1 8409 

.06908 

35.0 

-.26422 

-.77629 

3.75842 

1.86382 

1.05441 

.58408 

.29681 

.12897 

.04116 

40.0 

-.28802 

-.84314 

3.66096 

1.75287 

.94746 

.49500 

.23321 

.09156 

.02508 

45.0 

-.31080 

-.90653 

3.57030 

1 .65265 

.85477 

.42197 

.18479 

.06583 

.01562 

50.0 

-.33268 

-.96686 

3.48562 

1.56173 

.77403 

.36172 

.14760 

.04791 

.00994 

55.0 

-.35375 

-1.02441 

3.40631 

1 .47898 

.70341 

.31174 

.11884 

.03528 

.00646 

60.0 

-.37404 

- 1 .07937 

3.33194 

1.40350 

.64144 

.27009 

.09644 

.02630 

.00428 

M  =  10 


Ipo 

Uo 

a", 

a", 

a", 

a". 

all 

a\\ 

a" 

(in 

</:; 

10.0 

-.10578 

-.31631 

4.47621 

2.77375 

2.07820 

1.64130 

1.31465 

1 .04966 

.82610 

,63422 

.46880 

15.0 

-.13725 

-.40934 

4.32573 

2.57244 

1 .83486 

1.36648 

1.02012 

.74783 

.52944 

35464 

.21709 

20.0 

-.16438 

-.48887 

4.19922 

2.40748 

1 .64237 

1.15888 

.81033 

.54823 

.35089 

.20571 

.10349 

25.0 

-.18893 

-.56025 

4.08752 

2.26545 

1 .48228 

.99390 

.65301 

.40919 

.23778 

,12263 

.05082 

30.0 

-.21169 

-.62592 

3.98639 

2.13999 

1.34561 

.85918 

53165 

.30940 

.16382 

.07471 

.02570 

35.0 

-.23311 

-.68722 

3.89350 

2.02752 

1.22712 

.74738 

.43639 

.23642 

.11442 

.04638 

.01338 

40.0 

-.25343 

-.74493 

3.80738 

1.92573 

1.12336 

.65360 

.36074 

.18233 

08091 

.02929 

.007 1 7 

45.0 

-.27282 

-.79958 

.3.72708 

1 .8.3302 

1.03188 

.57435 

.30015 

.14184 

.05788 

.01881 

.00.395 

50.0 

-.29138 

-.85153 

3.65191 

1.74823 

.95082 

.50698 

.25128 

.11126 

.04188 

.01227 

.00223 

55.0 

-.30921 

-.90102 

3.58135 

1 .67043 

.87875 

.44945 

.21163 

.08799 

.03064 

.00814 

.001.30 

60.0 

-.32633 

-.94823 

3.51501 

1.59891 

.81448 

.40014 

.17931 

.07017 

.02268 

.00549 

.00077 

•S(cos  i^) 


V  (degrees) 


Figure  A.l  Stokes  Function  and  Modified  Kernel  According  to  Molodenskii  (^0  =  15°) 
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Figure  A.l  Stokes  Function  and  Modified  Kernel  According  to  Molodenskii  (y,„  15  ) 


RMS  Error  (meters) 


RMS  Error  (meters) 


kl 


\1  12 


Nl 

CAP  SIZE  ip„  (degrees) 


r 

i 


Figure  A. 3  Expected  Truncation  Error  (RMS)  for  Cieoid  Height  using  Method  III  and  Molodenskii's  Procedure. 
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Figure  A. 4  Expected  Truncation  Error  (RMS)  lor  Cieoid  Height  using  Method  III  and  Molodenskii's  Procedure  as 
Function  of  the  Number  of  Harmonics  Removed 
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Figure  A. 4  Expected  Truncation  Error  (RMS)  for  Geoid  Height  using  Method  III  and  Molodcnskii’s  Procedure  as  a 
Function  of  the  Number  of  Harmonics  Removed. 


Table  A. 2  Kxpected  Truncation  Krror  (meters)  for  (ieoid  Height  Based  on  Method  III  and  Molodenskii' 

Procedure 


d7v, 

5.0 

27.669661 

20.556343 

21.107032 

10.0 

21.739670 

14.294292 

15.528263 

15.0 

17.587212 

10.057129 

11.502553 

20.0 

14.777943 

7.337731 

8.555670 

25.0 

12.916628 

5.565207 

6.373636 

.30.0 

1 1.775998 

4.406253 

4.744508 

35.0 

11.229530 

3.798171 

3.521937 

40.0 

11.127562 

3.681 122 

2.602423 

45.0 

11.269953 

3.830107 

1.91 1084 

50.0 

11.471612 

4.043541 

1.392643 

55.0 

1 1.619539 

4.296507 

1.(8)5635 

60.0 

11.675333 

4.682291 

.718571 

65.0 

11.641289 

5.260989 

.507344 

70.0 

11.521156 

5.965933 

.35341 1 

75.0 

11.298.768 

6.642020 

.242489 

80.0 

10.9.79754 

7.131865 

.163589 

85.0 

10.416799 

7.330423 

.108290 

90.0 

9.728997 

7.204085 

.070178 

95.0 

8.917664 

6.789355 

.044406 

100.0 

8.066346 

6. 1 806(8) 

.027350 

M  4 

AN, 

<w,„ 

X! 

5.0 

27.669661 

7.143403 

8.1441  18 

10.0 

21.739670 

3.821151 

4.687022 

15.0 

17.587212 

3.916332 

2.780762 

20.0 

14.777943 

3.972130 

1.678309 

25.0 

12.916628 

3.379162 

1.021338 

30.0 

11.775998 

2.561030 

.622904 

35.0 

11.229530 

1.940179 

.379109 

40.0 

11.127562 

1.785319 

.229512 

45.0 

11.269953 

2.072378 

.137855 

50.0 

1 1 .471612 

2.514934 

.081967 

55.0 

11.619539 

2.845642 

.048144 

60.0 

1 1 .675333 

2.924490 

.027876 

65.0 

11.641289 

2.773443 

01*878 

70.0 

11.521156 

2.576939 

.008877 

75.0 

1 1.298368 

2.5563X7 

.004859 

80.0 

10.939754 

2.719752 

.(8)2*97 

85.0 

10.416799 

2.8*941  1 

(8)1  7*2 

90.0 

9.728997 

2.797542 

(H8)6X2 

95.0 

8.917664 

2.495368 

(88)717 

100.0 

8.066.346 

2.027741 

(8811  *6 
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Table  A. 2  (Continued) 


M  -  6 


d.\'( 

5.0 

27.669661 

3.264089 

4.202758 

10.0 

21.739670 

3.167549 

1.941453 

15.0 

17.587212 

3.645828 

.943511 

20.0 

14.777943 

2.9.36724 

.470071 

25.0 

12.916628 

1.892564 

.236672 

30.0 

1 1.775998 

1 .093629 

.119415 

15.0 

11.229530 

.640365 

.060057 

40.0 

11.127562 

.530456 

.029990 

45.0 

11.269953 

.685078 

.014823 

50.0 

11.471612 

.936809 

.00723 1 

55.0 

11.619539 

i. 308044 

.003473 

60.0 

11.675333 

1.698420 

.001638 

65.0 

11.641289 

1.900181 

.000756 

70.0 

11.521156 

1.826413 

.000341 

75.0 

11.298368 

1.622609 

.(XX)  149 

80.0 

10.939754 

1.551919 

(XXX)64 

85.0 

10.416799 

1.620535 

,(XXX)26 

90.0 

9.728997 

1.592219 

.(XXX)10 

95.0 

8.917664 

1.361978 

.(XXXXW 

100.0 

8.066346 

1.016394 

(XXXX)l 

M  -  8 

6N, 

5.0 

27.669661 

1.992536 

2.4^6855 

10.0 

21.739670 

3.037807 

924879 

15.0 

17.587212 

2.681047 

.371975 

20.0 

14.777943 

1 ,6<X)779 

.153939 

25.0 

12.916628 

.988105 

.064393 

30.0 

11.775998 

.747235 

.026968 

35.0 

11.229530 

.478502 

.011241 

40.0 

11.127562 

,38(XX)3 

(XU643 

45.0 

11.269953 

.570574 

.001  893 

50.0 

11.471612 

.794426 

.(XXI760 

55.0 

11.619539 

.883838 

(XX)299 

60.0 

11.675333 

.950892 

(XKH  l< 

65.0 

11.641289 

1.161459 

(XXKU3 

70.0 

11.521156 

1.309464 

OOOOI6 

75.0 

11.298368 

1  223234 

(KXKXK< 

80.0 

10.939754 

1,072946 

(XXXKI2 

85.0 

10.416799 

1  075179 

OOOOOI 

90.0 

9.728997 

1  .062328 

.(88X8X) 

95.0 

8.917664 

879693 

0000(8) 

l(X).() 

8.066346 

659448 

(XXXXX) 

.12 


Table  A. 2  (Continued) 


\/  10 


V"  (,v,  (,v,„  <*v„ 


5.0 

27.66966! 

1.750502 

1.542910 

10.0 

21.759670 

2.608021 

.477840 

1  5.0 

17.587212 

1.696188 

.159824 

20.0 

14.777945 

1.042415 

.055085 

25.0 

12.916628 

.944152 

.019181 

50.0 

11.775998 

.6151 20 

.006679 

55.0 

1 1.229550 

.296121 

.002511 

40.0 

11.127562 

.206595 

.(8)0791 

45.0 

11.269955 

.517089 

.(88)266 

50.0 

11.471612 

.486059 

.(888)88 

55.0 

11.619559 

.709787 

.(888)28 

60.0 

1 1.675555 

.798808 

.(8888)9 

65.0 

11.641289 

.784758 

.(8888)5 

70.0 

11.521156 

.899044 

.(8888)1 

75.0 

11.298568 

.949006 

.(88888) 

80.0 

10.959754 

.825826 

.(88888) 

85.0 

10.416799 

.780515 

.(88888) 

40.0 

9.728997 

.775971 

.(88888) 

95.0 

8.917664 

.625620 

.(88888) 

100.0 

8.066546 

.470106 

(88888) 

M  12 


4'" 

a- 

* 

. 

1 

«J<v„ 

5.0  ~~ 

27.669661 

1.716455 

1.015895 

10.0 

21.759670 

2.058058 

.260410 

15.0 

17.587212 

1.086172 

.072612 

20.0 

14.777945 

1.016599 

.020872 

25.0 

12.916628 

.757799 

.006057 

50.0 

1 1.775998 

.581679 

. (8)1756 

55.0 

11.229550 

.208457 

,(XX)505 

40.0 

1 1.127562 

.141440 

<XX)1 45 

45.0 

1 1 .269955 

.279572 

.(XXXMO 

50.0 

1  1.471612 

.410854 

.(XXX)  1  1 

55.0 

1 1 .619559 

.474925 

,(XXXX)5 

60.0 

1 1  675555 

.624421 

,(XXXX)1 

65.0 

II  641289 

.672848 

.(XXXXX) 

70.0 

1 1.521 156 

.650156 

.(XXXXX) 

7V0 

1 1.298568 

.728428 

(XXXXX) 

800 

I0.9597M 

.671656 

(XXXXX) 

85.0 

10.416799 

601482 

(XXXXX) 

900 

9  728997 

.596455 

(XXXXX) 

950 

891 "664 

.  4^44  5 1 

(XXXXX) 

100  0 

8.066546 

078855 

(XXXXX) 
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MODIFIED  VKNINCi  MEINESZ  KERNELS  AND  TRIM  ATION  ERRORS 


TABI.K  B.l  Coefficients  i+/„)  for  Modified  Veninj*  Meines/  Kernels 

Based  on  a  Molodenskii  Type  Procedure. 


M  2 


a\ 

a  i 

10.0 

.29863 

-4.50496 

15.0 

.47429 

-4.21927 

20.0 

.66262 

-3.91999 

25.0 

.85842 

-3.61841 

30.0 

1 .05679 

-3.32462 

35.0 

1.25385 

-3.04617 

40.0 

1.44691 

-2.78780 

45.0 

1 .63430 

-2.55183 

50.0 

1.81517 

-2.33875 

55.0 

1.98913 

-2.14784 

60.0 

2.15614 

-1.97770 

,Y/  =  4 

a\ 

a\ 

a\ 

10.0 

.29225 

-4.51548 

-2.82704 

-2.14386 

15.0 

.44868 

-4.26098 

-2.48351 

-1.72343 

20.0 

.60025 

-4.01985 

-2.16908 

-1.35669 

25.0 

.74279 

-3.79950 

-1.89448 

- 1 .05657 

30.0 

.87555 

-3.60103 

-1.66017 

-  .82044 

35.0 

.99943 

-3.42251 

-1.46178 

-  .63854 

40.0 

1.11574 

-3.261 17 

-1.2936S 

-  .49974 

45.0 

1.22569 

-3.1 1445 

-1.15052 

-  .39405 

50.0 

1.33022 

-2.98025 

-1.02798 

-  .31339 

55.0 

1 .42998 

-2.85693 

-  .92253 

-  .25155 

60.0 

1.52544 

-2.74326 

-  .83136 

-  .20.383 

M  6 

ci\ 

o' , 

a\ 

u\ 

a). 

10.0 

.28330 

-4.53025 

-2.84741 

-2.16952 

-1.74788 

1.43369 

15.0 

.41916 

-4.30910 

-2.54864 

-1.80345 

-1.32703 

-.97226 

20.0 

.54078 

-4.1 1517 

-2.29487 

-1.50595 

-1.004.3.3 

-.64394 

25.0 

.64954 

-3.94597 

2.08179 

1.26934 

-.76599 

-.42460 

30.0 

.74840 

-3.79614 

-1 .90075 

-1 .08000 

.59068 

-.28166 

35.0 

.83985 

3.661  17 

-1.74443 

.92631 

.46054 

.  1 8896 

40.0 

.92557 

3  S  3795 

1.60758 

.79984 

.36270 

-.12853 

45.0 

1 .00665 

3.42440 

1.48659 

.69462 

28827 

-  ,08874 

50.0 

1.08377 

3.31912 

-  1  ,37888 

.606.37 

.23108 

.06222 

55.0 

1.15737 

3.221 12 

1.28252 

.53190 

.  1 8675 

.04431 

60.0 

1.22776 

3.12968 

1.19606 

46876 

.15213 

.03205 

35 


Table  B.l  (Continued) 


M  =  8 

y>« 

lt\ 

u\ 

a\ 

a\ 

a\, 

a\ 

l/l 

10.0 

.27263 

-4.54786 

-2.87169 

-2.2001 3 

-1.78436 

-1.47550 

-1.22456 

-1.01102 

15.0 

.39001 

-4.35664 

-2.61304 

-1.88265 

-1.41855 

- 1 .07327 

-.79881 

-.57510 

20.0 

.49035 

-4.19615 

-2.40201 

-1.63354 

-1.14575 

-.79211 

-.52447 

-.32158 

25.0 

.57873 

-4.05762 

-2.22544 

-1.43396 

-.93949 

-.59504 

-.35039 

-.18090 

30.0 

.65896 

-3.93439 

-2.07326 

-1.26942 

-.77929 

-.45368 

-.23830 

-.10337 

35.0 

.73330 

-3.82248 

- 1 .9393 1 

-1. 13081 

-.65210 

.35002 

-.16465 

-.06024 

40.0 

.80307 

-3.71951 

-1.81979 

-1.01236 

-.54955 

-.27272 

-.11537 

-.03585 

45.0 

.86906 

-3.62400 

-1.71222 

-.91017 

-.46598 

-.21435 

-.08189 

-.02180 

50.0 

.93179 

-3.53494 

-1.61485 

-.82144 

-.39736 

-.16986 

-.05884 

-.01355 

55.0 

.99160 

-3.45161 

-1.52638 

-.74403 

-.34068 

-.13567 

-  .04279 

-.00860 

60.0 

1.04872 

-3.37348 

-1.44576 

-.67626 

-.29362 

-.10921 

-.03150 

-.00558 

M  =  10 

a\ 

u\ 

a\ 

a\ 

(I  l 

al 

a ' 

Ul 

al 

a  So 

10.0 

.26123 

-4.56667 

-2.89764 

-2.23283 

-1.82335 

-1.52021 

-1.27436 

- 1 .06520 

-.88273 

-.72160 

15.0 

.36354 

-4.39983 

-2.67160 

-1.95479 

-1.50205 

-1.16565 

-.89735 

-.67701 

-.49464 

-.34450 

20.0 

.44925 

-4.26227 

-2.48975 

-1.73849 

-1.26274 

-.91561 

-.64889 

.44165 

-.28246 

-.16348 

25.0 

.52470 

-4.14316 

-2.33625 

-1.56217 

-1.07638 

-.73189 

-.47925 

-.29526 

-.16551 

-.07889 

30.0 

.59337 

-4.03653 

-2.20222 

-1.41347 

-.92617 

-.59214 

-.35936 

-.20114 

-.09925 

-.03900 

35.0 

.65707 

-3.93919 

-2.08286 

-1.28546 

-.80249 

-.48344 

-.27260 

-.13908 

-.06071 

-.01981 

40.0 

.71685 

-3.84928 

-1.97526 

-1.17386 

-.69929 

-.39765 

-.20880 

-.09740 

-.03780 

-.01035 

45.0 

.77336 

-3.76561 

-1.87750 

-1.07577 

-.61239 

-.32927 

-.16134 

-.06902 

-.02.393 

-.00556 

50.0 

.82701 

-3.68738 

-1.78823 

-.98905 

-.53875 

-.27434 

-.12571 

.04946 

-.01540 

-.00370 

55.0 

.87811 

-3.61400 

-1.70641 

-.91207 

-.47603 

-.22995 

-.09875 

-.03585 

-.01007 

-.00174 

60.0 

.92686 

-3.54503 

-1.63125 

-.84352 

-.42240 

-.19389 

-.07821 

-.02628 

-,'X)669 
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Figure  B.l  Vening  Meinesz  Kernel  and  Modified  Kernel  Based  on  a  Molodenskii  Type  Procedure  (i+>„ 
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Figure  B.  1  Vening  Meinesz  Kernel  and  Modified  Kernel  Based  on  a  Molodenskii  Type  Procedure  (i^„  - 
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Figure  B.2  Vening  Meinesz  Kernel  and  Modified  Kernel  Based  on  a  Molodenskii  Type  Procedure  (if, 
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Figure  B.2  Vening  Meinesz  Kernel  and  Modified  Kernel  Based  on  a  Molodenskii  Type  Procedure  (if 


RMS  Error  (seconds) 


in)  K)  DEGREE  C  AP 


i  i  i  i  i  i  r  i 

0  2  4  A  X  III  i:  14 

NUMBER  OF  HARMONICS  REMOVED  (M) 


Figure  B.4  Expected  Truncation  Error  (RMS)  for  Vertical  Deflection  Using  Method  111  and  a  Molodenskii  Type 
Procedure  as  a  Function  of  the  Number  of  Harmonics  Removed 
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Figure  B.4  Expected  Truncation  Error  (RMS)  for  Vertical  Deflection  Using  Method  II’  and  a  Molodenskii  Type 
Procedure  as  a  Function  of  the  Number  of  Harmonics  Removed, 


TABLE  B.2  Expected  Truncation  Error  (seconds)  for  Verlical  Defieclion 
Based  on  Method  III  and  a  Molodenskii  Type  Procedure 


M  -  2 


CAP  SIZE 

1^0 

<50, 

5.0 

3.978880 

3.662607 

3.663817 

10.0 

3.111038 

2.742907 

2.752567 

15.0 

2.505130 

2.096936 

2.126250 

20.0 

2.031197 

1.590299 

1 .648072 

25.0 

1.652192 

1.186509 

1.273362 

30.0 

1.352733 

.873341 

.9788(8) 

35.0 

1.121781 

.644263 

.748134 

40.0 

.947027 

.489385 

.568485 

45.0 

.814783 

.390303 

.429320 

50.0 

.712961 

.323158 

.322058 

55.0 

.634082 

.269139 

.239799 

60.0 

.575874 

.221037 

.177051 

65.0 

.539023 

.182188 

.129476 

70.0 

.523560 

.161844 

.093660 

75.0 

.526517 

.166677 

.066921 

80.0 

.542567 

.193653 

.047152 

85.0 

.566245 

.235818 

.032702 

90.0 

.593267 

.288449 

.022278 

95.0 

.620287 

.347719 

.014873 

100.0 

.643988 

.408531 

.009703 

M  =  4 

CAP  SIZE 

<50, 

<50  „ 

5.0 

3.978880 

2.490406 

2.496183 

10.0 

3.111038 

1.469157 

1.508087 

15.0 

2.505130 

.842184 

.931024 

20.0 

2.031197 

.479325 

.575539 

25.0 

1.652192 

.364583 

.355742 

30.0 

1.352733 

.379008 

.219870 

35.0 

1.121781 

.383924 

.  1 35749 

40.0 

.947027 

.347850 

.083569 

45.0 

.814783 

.283567 

.051 187 

50.0 

.712961 

.211813 

.031125 

55.0 

.634082 

.148805 

.018749 

60.0 

.575874 

.101453 

.01 1 165 

65.0 

.539023 

.069449 

.006559 

70.0 

.523560 

.053782 

,003792 

75.0 

.526517 

.056723 

.002153 

80.0 

.542567 

.070400 

.001 197 

85.0 

.566245 

.088669 

.(88)650 

90.0 

.59.3267 

.113041 

.(88)344 

95.0 

.620287 

.143649 

(88)176 

100.0 

.643988 

.174213 

.(888)87 

Je. 
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Table  B.2  (Continued) 


6 


CAP  SIZE 

66, 

<50,, 

5.0 

3.978880 

1.855658 

1.869405 

10.0 

3.111038 

.841762 

.913194 

15.0 

2.505130 

.377418 

.454259 

20.0 

2.031197 

.321540 

.227998 

25.0 

1.652192 

.343545 

.115407 

30.0 

1.352733 

.299329 

.058705 

35.0 

1.121781 

.218555 

.029874 

40.0 

.947027 

.154355 

.015147 

45.0 

.814783 

.131300 

.007625 

50.0 

.712961 

.120768 

.003799 

55.0 

.634082 

.099153 

.001869 

60.0 

.575874 

.069989 

.000905 

65.0 

.539023 

.044329 

,(XK)430 

70.0 

.523560 

.030527 

.000200 

75.0 

.526517 

.033051 

.(XXXWl 

80.0 

.542567 

.043041 

,(XXX)40 

85.0 

.566245 

.052876 

.(XXX)  17 

90.0 

.593267 

.066097 

,(XXXX)7 

95.0 

.620287 

.085426 

,(XXXX)3 

100.0 

.643988 

.102787 

.(XXXX)l 

=  8 


CAP  SIZE 

<50, 

<50,,, 

<56,, 

5.0 

3.978880 

1.419918 

|  .444878 

10.0 

3.111038 

.484336 

.571066 

15.0 

2.505130 

.283755 

.231284 

20.0 

2.031197 

.304085 

.095652 

25.0 

1.652192 

.238663 

.040132 

30.0 

1.352733 

. I 57320 

.016938 

35.0 

1.121781 

.138061 

.007142 

40.0 

.947027 

.135345 

,(X)2994 

45.0 

.814783 

.107694 

.(X)  1 243 

50.0 

.712961 

.073710 

(XX)509 

55.0 

.634082 

.055312 

.(XX)205 

60.0 

.575874 

.043586 

,(XXX)8 1 

65.0 

.539023 

.028934 

,(XXX)3 1 

70.0 

.523560 

.018118 

.(XXX)  12 

75.0 

.526517 

.020591 

.(XXXXM 

80.0 

.542567 

.029886 

,(XXXX)2 

85.0 

.566245 

.037026 

.(XXXX)I 

90.0 

.593267 

.045 1 38 

(XXXXX) 

95.0 

.620287 

.058739 

(XXXXX) 

100.0 

.643988 

.068533 

.(XXXXX) 

Table  B.2  (Continued) 


M  10 


'AP  SIZl: 

69, 

«*8W 

5.0 

3.978880 

1.094937 

1 . 1 33355 

10.0 

3.111038 

.303695 

.363392 

15.0 

2.505130 

.279978 

.121074 

20.0 

2.031197 

.226822 

.041608 

25.0 

1.652192 

.140284 

.014535 

30.0 

1.352733 

.133846 

.005102 

35.0 

1.121781 

.122272 

.001786 

40.0 

.947027 

.0861 17 

.000620 

45.0 

.814783 

.071222 

.00021.3 

50.0 

.712961 

.063705 

.000072 

55.0 

.634082 

.044521 

.(XX3024 

60.0 

.575874 

.028111 

.000008 

65.0 

.539023 

.018260 

,(XXXX)2 

70.0 

.523560 

.010516 

.(XXXX)l 

75.0 

.526517 

.012848 

.OtXXXX) 

80.0 

.542567 

.021672 

.(XXXXX) 

85.0 

.566245 

.028007 

.000000 

90.0 

.593267 

.033488 

.(XXXXX) 

95.0 

.620287 

.04366.3 

.000000 

100.0 

.643988 

.048752 

.(XXXXX) 

M  =  12 


ap  si/.i; 

4>„ 

69, 

5.0 

3.978880 

.843845 

.896296 

10.0 

3.111038 

. 248 1 06 

.234207 

15.0 

2.505130 

.247209 

.064694 

20.0 

2.031197 

.145935 

.018544 

25.0 

1.652192 

.125036 

,(X)5402 

30.0 

1.352733 

.113842 

.001579 

35.0 

1. 121781 

.078850 

.000459 

40,0 

.947027 

.074463 

(XK)1 32 

45.0 

.81478.3 

.060512 

,(XXX)37 

50.0 

.712961 

.042577 

.(XXX)  10 

55.0 

.634082 

.035972 

,(XXXX)3 

60.0 

.575874 

.024122 

.(XXXX)l 

65.0 

.539023 

.012894 

.(XXXXX) 

70.0 

.523560 

.006329 

.(XXXXX) 

75.0 

.526517 

.008318 

.(XXXXX) 

80.0 

.542567 

.016099 

.(XXXXX) 

85.0 

.566245 

.022160 

(XXXXX) 

90.0 

.593267 

.026173 

(XXXXX) 

95.0 

.620287 

.034066 

.(XXXXX) 

100.0 

.643988 

.036357 

(XXXXX) 
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APPENDIX  C 


DERIVATION  OF  CERTAIN  INTEGRALS  INVOLVING  LEGENDRE  FUNCTIONS 


DERIVATION  OF  CERTAIN  INTEGRALS  INVOLVING  I.KGKNDRK  FUNCTIONS 


Use  will  he  made  of  the  reeurrenee  relations  as  given  hy  I  ebedes  (1972).  In  all  eases  the  integral  ion  inter¬ 
nal  will  he  -  I  <  u  <  v  <  b  <  I. 

I.  Definitions: 

P:(.v)  =  (l  v’)-  <r  P\ ,v) 
dv 

PAx)  =  d  PAx)  (CM) 

(lx 


II.  Reeurrenee  relations: 

( i  -  ,v* )  d  pzz  =  («  +  »i )  r:  i  -  nx  r: 

(lx 

nP„=xP„  -  P.  , 
hP„-,  =  P„  -xP'„.\ 

III.  Theorems: 

(i)  Theorem  C-l  : 

(n  -  A)  xP:p7  -  (n  +  m)  PZ  , PZ  +  (A'  +  hi)  PZPZ  , 
(n  -  A )(//  +  A  +  I) 

for  «  *  A  and  n  >  I ,  A  >  1 . 


f 


/j:(a)  r/.v  = 


(C.2) 

(CM) 

(C.4) 


(C.5) 


Proof:  Consider  the  following  differential  equations  satisfied  hy  the  associated  Legendre  functions 

d  {(I -vVt  +  {n(n  t  I)-  "'‘.^  =  0  ((h) 

r/.v  I  -  vJ 

with  a  solution  given  by  PZ(x),  and 

d  f(l  -,v')v'»  +  {A(A  t  I)-  m‘ }  v  =  0  (C.7) 

(lx  I  -  V! 

with  one  of  its  solutions  PZ(x).  Assume  that  n  *  A  and  multiply  equation  (C.6)  hy  P?(x)  and  (C.7)  hy  PZ(x ), 
form  their  difference,  and  then  integrate  to  get 


r  pr/r/  (ix 

■')  r/  *»:] 

.  r/.v  -  J 

f  p:  d  - 

v'»  d  p: 

•  \(lx 

r/.v  J 

r/.v  1 

(lx 

*  \n(n  t  I)  A(A  ♦  I )\  J  PZPZdx  0.  (C.K) 


) 


* 
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Integrating  the  first  two  integrals  by  parts  and  noting  that 
//(/;  +  1)  -  A(A  +  I )  =  (n  -  k)(n  +  k  +  I) 
leads  to 


f  P?P7fix=  0-x!) 

jp;  d  p:  -  p:  d  P7 

L  (n  -  A)(n  +  A  +  1) 

1  dx  dx 

(C.9) 


Using  (C.2)  and  rearranging  terms,  gives 


p:  dx  = 


in  -  k)  xP"PT  -  in  +  m)  P"-,P?  +  (A'  +  m)  P, 


in  -  k)(n  +  k  +  1) 


for  n  *  k  and  n  >  ! ,  A  >  1 . 


(ii)  Theorem  C-2: 


f  Pi  Pi  dx  =  -  nk 

in -k) xP,P„  -  in  +  1)P„P,.,  +  (A  +  l)P„.,Pt  ' 

* a 

in  -  A)(«  +  A  +  1) 

for  n  *  A  and  n  >  2,  A  >  2. 

Proof:  Set  in  =  I  in  equation  (C.5)  to  get 

in  -  k)xP\Pl  -  in  +  1)  Pi,  P\  +  (A  +  I)  PlPj, 
(«  -  A)(n  +  A  +  I) 

From  equation  (C.l)  for  m  =  I,  and  (C.3) 
p]  =  (i  -x‘y'‘p; 

Pl,=  (I  -  .V1)1"  P’., 

=  (I  -  .v')1'1  (.vP,  -  jP ,),  for./  S*  2. 


s. 


Pi  Pi  dx 


Substituting  the  above  into  (C.ll)  and  simplifying  yields 


* 


/  n(n  f  1 )  PJ\  -  A(A  +  I  )  P„P. 
\  in  -  A  )in  f  k  +  I ) 


(C.10) 


(C.ll) 


((12) 
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Finally,  using  equation  (C.2)  with  in  -  0,  gises 


/  Pi.  Pi  dx  =  -  nk 

(n  -  k)  xP„Pt  -  (n  -t  \)  PJ\  ,  i  (k  t  1 )  P„.,Pi 

* 

(m  -  A  )(M  +  A  4  1 ) 

tor  n  *  A  and  n  >2  ,  k  >  2. 


(iii)  Theorem  C-3: 


f  P.(.v)  P,(.v )  r/.v 

(m  -  A )  xP„P,  -  nP..  ,Pt  +  A  P„P,  , 

L  (m  -  A)(m  +  A  +  1) 

for  n  ±  k,  n  >  I ,  k  >  1 . 

Proof:  Set  in  =  0  in  equation  (C.5)  and  the  result  follows. 

(iv)  Theorem  C-4: 


f  Pl  dx  =  f  PL.  dx  + 

A(P'  +  P'.,)  -  2P..P..., 

(2m  +  1) 

(2m  +  1) 

Proof:  Replacing  one  of  the  /\,(.v)’ s  by  (C.3)  gives 
nf  P‘cix  =  J  xP„P„  dx  -  J  P„P fix. 

Integrating  by  parts  the  integrals  on  the  right  side  of  (C.15)  give 

J.  xP.P.dx  =  l  [.vT-i:  -  I  j  Pi  dx 
and 

J  P.P..,dx  =  [/>„/>„.,  l,-  J  P»  i P«  dx 

or  using  identity  (C.4) 

J  PJ\.,  dx  \PJ\.,\"  -  n  J  Pi  ,  dx  -  J  xP,  ,/M,  dx. 

Substituting  (C  .  16)  and  (C.I7)  into  (C  .15)  and  simplifying  gives 

(2«  *-  I )  f  Pldx  \xP‘.-lPJ\X  *  2m  J  Pl,dx  +  2  /  \P„  ,/\,  ,  dx. 
Integrating  the  last  integral  of  (C.I8)  by  parts  gives 

fmxP.,P.,dx  l\xPl,\[  [  J  Pi ,, 


(CM3) 


(C  . 14) 


(C.15) 


(C.  16) 


dx 


(C  .  17) 


(CIS) 


(('.111 
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(C.20) 


so  equation  (C.18)  can  be  written  as 

f1'  b 

(2n  +  1)  J  P‘dx  =  (2 n  -  1)  J  Pi,  dx  +  [x(Pl  +  Pi,) -2P.P..,]  t 


or 


f  Pl  dx  =  <2”-  1}.  /  Pi,  dx  + 
(2  n  +  1)  ° 


x(Pl  +  Pi,)  -2 PnP, 


(2  n  +  1) 


(v)  Theorem  C.5: 

rh  rh 

J  [Pi]2  dx  =  n(n  +  1)  J  Pi  dx  +  n[P„P„.,  -  xPl\\  forn»2. 
Proof:  From  (C.l)  with  m  =  1  and  (C.2)  with  m  =  0 
[P^]2  =  Pl(l  -*2)/\ 

=  P'„(nP„~,  -  nxPA. 


f 

^  d 

Using  (C.4)  and  then  (C.19)  in  the  first  integral  on  the  right  of  (C.22)  yields 

Jo  P*-,P 

=  (n-bi  Pl,dx+±[xPtX. 

Using  equation  (C.14)  in  equation  (C.23)  gives 
J  P.  ,Pldx  =  (n  +  {)  Jo  Pldx+  [/>„/>„-,  -  i xPll  . 

Finally,  substituting  equations  (C.I6)  and  (C.24)  into  (C.22)  and  simplifying  yields 
J  [Pi]1  dx  =  n(n  +  1)  J  Pi  dx  +  n\PrP„.,  -  xPt )„  . 


,  dx 


n  J  Pi.  dx  +  xP„.,P„ 


dx 


[Pi] 2  dx  =  n  r  P.-,P.  dx  -  n  xP*P„  dx. 
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a* 


(C.21) 


(C.22) 


(C.23) 


(C.24) 
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